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SUMMARY 


Inconsistencies between analytical results and measurements for 
liquid rocket thrust chamber performance, which escape suitable 
explanations, have motivated the examination of the basic 
physical modeling formulations as to their unlimited application. 
The publication of Prof. D. Straub's book, "Thermofluid-dynamics 
of Optimized Rocket Propulsions", further stimulated the interest 
of understanding the gas dynamic relationships in chemically 
reacting mixtures. A review of other concepts proposed by Falk— 
Ruppel (Gibbsian Thermodynamics), Straub (Alternative Theory, 

AT) , Prigogine (Non-Equilibrium Thermodynamics) , Boltzmann 
(Kinetic Theory) , and Truesdell (Rational Mechanism) has been 
made to obtain a better understanding of the Navier-Stokes 
equation, which is now used extensively for chemically reacting 
flow treatment in combustion chambers. 

In addition to the study of the different concepts, two workshops 
were conducted to clarify some of the issues. The first workshop 
centered on Falk-Ruppel ' s new "dynamics" concept, while the 
second one concentrated on Straub's AT. 

In this report brief summaries of the reviewed philosophies are 
presented and compared with the classical Navier-Stokes 
formulation in a tabular arrangement. Also the highlights of 
both workshops are addressed. yr - 

££-ri :A 

The primary conclusions are as follows: 

!• The new thermodynamics theory by Faulk— Ruppel provides 
an indepth approach to formulate thermo-fluid equations. 

However, the connection with the conventional Navier-Stokes 
equation is not apparent. Further investigation of this issue is 
necessary for the assessment of the governing equations in 
combustion chamber flows for CFD simulations. 

2. Both, NASA's ODE program and the Munich Method Code for 
thermodynamic equilibrium calculations of rocket engines now 
carry an option to consider finite area combustion effects. The 
solution process in both approaches is based on the Lagrangian 
m ^^^"4.plier scheme. The results from these two programs are 
practically identical for two demonstrated sample cases. 

However, the ODE program contains many more chemical elements and 
species and is equipped with several practical calculation 
options. 

3. The derivation of a new "Navier-Saint Venant" equation 
as a result of Prof. Straub's "Alternative Theory" was not 
P**®^®Dted. Hence, the equation was not accepted as a preferred 
scheme over the Navier-Stokes Equation at this time. Some 
interesting results should be surveyed. 

4. The comparison of Navier-Stokes equation with other 
formulations from Boltzmann, Prigogine, and Truesdell have 
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revealed certain limitations, when this equation is used for 
combustion chamber flow modeling. 

5. This review process is an initial attempt to establish a 
sound foundation for the analytical simulation of rocket engine 
flow processes. Further re-evaluations are in order and should 
be supported by a group of experts from the academic profession 
to assist the engineering community in solving their problems 
with a high degree of confidence. 
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I . INTRODUCTION 


Since the introduction of the maximization of entropy for 
equilibrium in rocket engine combustion processes by Prozan [1] 
in 1969, this scheme has induced considerable different opinions 
in the JANNAF community. In 1982, NASA/MS FC awarded a contract 
(NAS8-34946) to Continuum Inc. to program this entropy 
maximization concept and provide quantitative evidence. The 
results led to a workshop at Continuum Inc. , where Prozan / s 
proposed concept was reviewed and debated. In the meantime, 
some new developments appeared in Europe. The introduction of a 
"Dynamics" concept by Falk and Ruppel [2,3] forms the basis for a 
new look at the currently used thermodynamic and fluid dynamic 
laws. Building on their concept, Prof. Straub formulated an 
"Alternative Theory (AT) " which is characterized in a recently 
published book entitled "Thermofluid-Dynamics of Optimized Rocket 
Propulsions" . The contents of the book became strongly 
controversial. To gain further understanding of these 
differences, MSFC awarded a contract to The University of Alabama 
m Huntsville (UAH) to investigate the fundamental analytical 
formulations, which are currently used to simulate the physical 
processes in liquid rocket combustion chambers and nozzles. 

To accomplish these goals, two workshops were organized, and 
their proceedings have been published [4,5]. Brief information 
about both workshops and the conclusions are included in Appendix 
I and II of this report. The foundation of the "Gibbsian 
Thermodynamics" and one important application’ is described in 
Section II. Information about Prof. Straub's (AT) is included in 
Section III. Section IV discusses the derivation of conservation 
laws via Boltzmann equations, and the ones based on Prigogine’s 
non-equilibrium thermodynamics are presented in Section V. 
Truesdell's rational mechanics theory is displayed in Section 
VI. a comparison of these concepts and theories with concluding 
remarks are presented in Section VII. 


The motivation to select the classical theorems (Gibbsian, 
Boltzmann, Prigogine and Truesdell^- to be included in this study 
is because in one way or another they attempt to formulate the 
non-equilibrium chemically and radiatively acting flow. These 
are being compared with the current Navier Stokes Equation used 
in the present Computational Fluid Dynamics (CFD) for the 
combustion chamber flow simulation. 
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II. GIBBSIAN THERMODYNAMICS 


Professors Gottfried Falk and Wolfgang Ruppel are the originators 
of a new concept, called "Gibbsian Thermodynamics", which is 
presented in the following two books [2,3], written in German and 
printed by the Springer Publishing Company: 

1. Mechnik, Relativitat, Gravitation - (1977) 

2. Energie und Entropie - (1976) 

Some important chapters of these books have been summarized and 
translated, and they are included as Appendix C in [4], This 
section provides only excerpts from [4], relating to the 
justification of the new concept and one notable application. A 
comprehensive document showing the derivation of the novel 
equation and its advantage should be prepared in English to 
assure proper attention and relevant studies. 

In the book "Energy and Entropy" the mass character of 
energy, the various forms in which energy is exchanged, energy 
fluxes, and the spatial distribution of energy or energy flux 
concentration are addressed. The mass character of energy 
permits the summation of all energy forms. The energy which is 
connected with an energy transport during a process appears in 
specific energy forms. Please note that only the transported 
energy relates to an energy form and not the stored energy. 

Energy forms are: energy of rotation, energy of motion, energy 

of compression, surface energy, chemical energy, heat energy, and 
others. All energy forms can be represented with the following 
mathematical description f dX, where X is the extensive variable 
specifying the energy form and f is the intensive variable, 
identifying the quantity of energy which is exchanged in this 
form. 

Energy Form = (intensive variable f) • d (extensive variable X) 

There are as many different energy -forms as there are extensive 
variables, however, they do not have to be independent of each 
other. The number of independent energy forms reflect the number 
of degrees of freedom in a system. The general validity of this 
energy form convention cannot be proven, like no assertion about 
nature can be confirmed. 

One of the great achievements in physics has been the 
identification of "standard variables", which permit the 
description of an untold number of processes in nature. A 
special simple treatment of such processes can be accomplished 
with variables which satisfy conservation laws. They are 
especially favored as standard variables. 
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Since energy is a mass guantity, which is distributed in 
space and can flow, the following Gibbs Fundamenta I Form is 
introduced. The terminology "form" is used because the~~equation 
is composed of differentials of the first order (linear) also 
known as Pfaff's form. 


dE - i” i dX j + f 2 dX 2 + . .. + £ n dX n (II. 1) 

Distinct energy forms are presented below with bold letters 
identifying vector quantities: 

Energy of : 


ROTATION: 

ftdL 

TRANSLATION: 

-Fdr 

MOTION: 

vdP 

COMPRESSION: 

-pdV 

SURFACE : 

crdA 

ELECTRICAL: 

pdQ 

CHEMICAL: 

udN 

HEAT: 

TdS 

POLARIZATION: 

Ed 7 

MAGNETISM: 

Hdm 


(Angular Velocity) * L (Rot. Momentum) 
(-Force) * r (Position Vector) 

(Velocity) * P (Momentum) 

(-Pressure) * V (Volume) 

(Surface Tension) * A (Surface Area) 
(Electrical Potential) * Q (Charge) 
(Chemical Potential) * N (Partical No.) 
(Temperature) * S (Entropy) 

(El. Field Strength) * 7 (Dipole Moment) 
(Mag. Field Strength) * m (Dipole Moment) 


Gibbs Fundamenta] Form 


dE-Ef n dX n (II. 2) 

In this equation dE is the received or produced energy. It 
is important that each energy form is independent of the other 
ones, which means that the energy flux of each energy form must 
be possible regardless of the other fluxes. 


If all extensive variables have a constant value, then all 
differentials are zero, which in turn requires dE = 0 or E = 
constant. The corresponding Gibbs Function reads: 

E = f (P, r, L,S,V,A,N) (II. 3) 

(The electrical, polarization, and magnetic energy forms are not 
considered here) 
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This function describes systems which are physically equal 
when the same extensive variables are involved. However, they do 
not need to be the same in reality. 

Gibbs functions are also called Thermodynamic Potentials. 

Statement 


The identification of the energy E of a system as a function 
of the independent extensive variables is a very difficult task, 
which can only be solved by approximations in the end. Large 
disciplines of physics, such as statistical mechanics, are used 
primarily to find these functions. 

The intensive variables are obtained by differentiating the 

Gibbs function E = E ( ) with respect to the extensive 

variables . 

3E 3E 3E 3E 3E 

dE = — dP - — dr + — dS + — dV + — dN (II. 4) 

3P 3r 3S 3V 3N 

where the partial derivatives represent the velocity v = 3E/3P, 

the force F = 3E/3r, the temperature T = 3E/3S, the pressure -p = 
3E/3V, and the chemical potential u = 3E/3N respectively. These 
results give 


dE = vdP + Fdr + TdS + pdV + udN (II. 5) 

The specific features of this equation are: 

o The fundamental equation reflects an energy change and not 
its conservation. 

o All energy terms are portrayed by the same mathematical 
structure. 

o All terms have quantity character and become additive. 

o Each term represents a specific energy form. 

o Each energy form is presented by two conjugate physical 
parameters of an intensive and extensive variable. 

o Each variable can appear in the equation only once. 

o A system can be subdivided if the variables in one subsystem 
do not affect the variables in the other one. 

The identified "Gibbs Fundamental Form" and its identified 
features form the basis of Gibbsian Thermodynamics. In order to 
show the connection with the classical mechanics (i.e. Newton's 
Second Law) we shall discuss the derivation of Newton's Second 
Law from the Gibbs Fundamental Equation. Recall the Gibbs 
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fundamental equation. 


dE = i dX j 

E = E(L,r,P,V,A,Q,Nj ,S,7,m) (II. 6) 

Separating the energy into two parts which represent the motion 
and the internal character, neglecting the energy forms of 
polarization and magnetism results in: 

E(P,r,L,S,V,A,N; ) = E : (P,r) + E 0 (L, S , V, A, N j ) (II. 7) 

Reducing this equation to hold only the energy of motion (vdP) 
and energy of translation (Fdr) requires that the dynamic 
variables of the internal energy subsystem have: (1) constant 

values and, therefore, their differentials are zero, (2) their 
differentials are not zero, but the summation of the associated 
energy forms are zero, (3) the independent variables 
representing the internal energy have no influence, or perhaps 
only a weak one on the other subsystem. 


i^dX! = vdP; Velocity • d (Impulse) 

f 2 dX 2 = Fdr; -Force • d (Position) (II. 8) 

With a further assumption that the system does not exchange 
energy with another system finally results in the equation: 


dE = vdP - Fdr = 0 

or (II. 9) 

E = E(P,r) = const. 

In a static field where the force is only a function of the 
position, F = F(r), which postulates that the velocity is only a 
function of the impulse, v = v(P) , the differentials of the 
extensive variables become total differentials. The energy of 
motion is now representing the kinetic energy, and the energy of 
translation the potential energy. — 

If the system is subject to a change in time and postulating 
that P and r are functions of time, then differentiation with 
respect to time yields: 


dP dr 

v — - F — =0 (II. 10) 

dt dt 

Recognizing that the change of position with time is equivalent 
to the kinematic velocity, v = dr/dt, which can be equated for 
point particle motion with the dynamic velocity v = 3E/9P, and 
substituting the proper term in the equation leads to: 
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dP 

v — -Fv = 0 
dt 

(II. 11) 

'dP 

v F =0 

\dt 

This equation is satisfied when 

v = 0 

'dP 

F =0 

vdt ) 

and the scalar vector product is Zero. 

The second requirement leads to the equation below 

dP 

— = F (11.12) 

dt 

Please note that there are two equations of motion which are 
valid in general when only the energy of motion and energy of 
translation are present, where F = F(P,r) and v = v(P,r). 
Separation into two independent subsystems requires that they 
cannot share independent variables, and total derivatives can be 
used. 

dP 

— = F 
dt 

(11.13) 

dr 

dt 

Now, the intensive variables F and v are only functions of either 
P and r respectively. 

For Newton's condition (v « c) the two equations translate 
into a better known form. With the assumption that a field of 
the first type is present, the force F is only a function of r, 
which requires that the velocity v is only a function of P. The 
latter condition can be obtained in the following approach. In 
Einstein's important relationship the impulse is related to the 
energy in the following way: 

E 

p = — v (11.14) 

c 4 
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In this equation E represents the total transport energy which 
itself depends on v and increases with growing velocity. For 
transport velocities v 4 0 the energy approaches a finite value 
which is known as the internal energy E 0 = M c 1 2 . Substituting 
this term in the previous equation results in 


This is just the relationship we are looking for, where v is a 
function of P. In this equation M represents the inert mass, 
which is normally a function of the velocity, but here M is 
equivalent to the velocity independent rest energy or internal 
energy E 0 for v = 0 and P = 0. 


Substituting the impulse value for Newton's condition (v « c) in 
the bracket term (II. 11), and taking M out of the differential 
expression, since it is a constant and therefore independent of 
time, yields the following equation: 


dP 

— = F 
dt 

d (Mv) 

= F 

dt 

dv 

M — = F 
dt 

d 2 r 

M = F 

dt 2 


(11.16) 


The last equation is Newton's famous law: 


MASS * ACCELERATION = FORCE 


During this derivation the following assumptions have been made: 

o The motion of the system is virtually independent of the 
internal energy of the system with v << c or c|p| << E. 

o The energy of motion and translation of a system remain 
when : 


(1) The fixed values of the extensive variables L 

(Rotational Momentum) , S (Entropy) , V (Volume) , N (Particle 
Numbers) and therefore, also their conjugated intensive 
variables A (Angular Velocity) , T (Temperature) , p 
(Pressure), u (Chemical Potential), state that the internal 
energy does not change. This relates to the motion of rigid 
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bodies or mass particles. 

(2) The sum of the energy forms representing the internal 
energy is zero, but the individual terms are not. Here the 
motion is independent of the internal energy. 

f 3) The independent variables portraying the internal 
energy do not affect the other subsystem or have only a weak 

impact on it. 

The mass M is independent of time. 

Only two forms of energy are considered which are 
independent of each other. 

The transitional energy Fdr is restricted even to a static 
field, where F = F(r) and not of v. 


Please note: 
fields) with 
relationship 
cannot result 


If fields of the second type (non-conservative 
F = F(P,r) and v = v(P,r) are present, the 
p = T«(v is not possible. Therefore, equation (11.16) 
for such a condition. 
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III. ALTERNATIVE THEORY: STRAUB METHOD* 

Straub has deduced his so-called "Alternative Theory" (AT) , 
[6], on the basis of the 'dynamics' theory by Falk and Ruppel, 
[2,3], combined with Progogine's new micro-theory [7], The 
Navier-Saint Venant equation is a principle result of the AT and 
leads to a confrontation with the Navier-Stokes equation of 
motion. Subsequently, Straub's description of the AT from [6] is 
presented. 


III.l. The Navier-Saint Venant Equation of Motion 

In order to offer a quantitative presentation of the 
numerous important results of the AT, the equation of motion of a 
compressible multicomponent single-phase fluid mixture is 
compared to the Navier-Stokes equation of motion. The basic 
relationship representing the momentum conservation law is a 
Cauchy-type equation of motion. 

pDV = pf - 7 • 7T + 1/2 9 tP *, (III.l) 

whereas D is a differential operator called the substantial or 
material time derivative 


D: = 9 t + v • 7 (III. 2) 

in this definition 9 t denotes the partial time derivatives and 7 
is the gradient vector operator (nabla-operator) . 

The local flow velocity v involved has a double definition: 
the original 'dynamic' establishment 


9E(P,r,s,V,N k ) 

v = 

9P 


(III. 3) 


is obtained with the Gibbs Fundamental Equation [6] of the 
multicomponent single-phase fluid mixture. Identification of v 
with the kinematic velocity v*: = dr/dt guarantees conservation 
of the particle number of all baryons. This invariance 
generalizes the Newtonian mechanics of finite systems of mass- 
points, called 'bodies' . At every given instant the body is 
located at some place and assigned real numbers denoting the 
properties of state. In field theories the continuity equation 
controls the time behavior of such a body by considering the 
baryons constancy (see Straub 1988, pp. 53 f) . 


*The content of this section is directly transposed form ref. [61. 
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with^mass^density ! and l appearin 9 in equation (in „ 

gravity), are defied fs foii SPeCi “ c fie ^ force J fioi a " g 

ea as follows in the AT: ff r exaiT1 Ple, 


7r ‘ P5 1 Pressure tensor 
t: = v - i dissipation velocity 


(HI. 4) 
(HI. 5) 


Here , stands for the unitv t " (11 —' 

o? e «st e r If ab TMs e ^ With theS ^atiSn ?or S the e h therm ° dynan,ic 

Wlth the th -“i eqi S r oV; t a P tT ltted to be 


p * R T p Z 


where p denotes the density j (III. 6] 

? hypothetical state ot ?est ih the tei "P*rature defined f 
P) 13 ldenti = a l to one only^it^i^flfliWlity facto? z l 
of viscous stress tensor +. 

di “ i classes 

T ' 1/2 pV * viscous stress tensor 
The quantity * equals rh (HI. 7) 

Vel ° City * - --- -cai 

. v ^e i, the limiting 

lim <f> = o 

cr-> 0 

disappears only for va nj v. • (III. 8) 

values^of^ the^ entropy by the iocal 

to t h e F °c r iLs = i?ir £au°’ h j equation*?? <&£ ^™aUy corresponds 

° f m ° tion '-r a 

F - *“ st ress P tensor ° t 


12 


s l 1 I 

r = t T Z" 1 p - < D - — (7 • V) 8 - | 

R \ 2 j I 

! 

1 ! 

t T T ( 3 1 s + 7 < j s ) 8 - i 

2 i 

i ,1 (III - 9) 

- t T ( ZR) ~ 1 p] 7SV - - (7 • sv) 8 | 

\ 2 j 


- t T (ZR ) -1 p] [vs xv] x s| 


is theoretically derived. 

For simplification, the abbreviation 


s 

t T pTs : = 25 = t T Z -1 p - (III. 10) 

R 

is introduced. 

Before equation (III. 9) is discussed in greater detail, its 
corresponding traditional continuum mechanics relationship should 
be offered for comparison. These Navier-Stokes equations, as 
generalised notations of the set of field equations for v (r,t), 

i P( r it) and p(r,t) , are defined by its viscous stress 
tensor: 


t: = 2uD + i) v - — ttj 8 7 • v (III. 11) 


Perhaps the clearest interpretation of this equation is found in 
Truesdell (1984, pp. 426-427 & 409-410). In this case, D is the 
deviator which is the symmetrical share of the (tensorial) 
velocity gradient 7V in symmetrical stress tensors; u and ^ are 
the fluid's shear or volume viscosity respectively. c. 

Truesdell offers convincing reasons why the known Stokes relation 
(i.e. y) v e 0) is valid. 

At first glance one notes that the expression (III. 9) for 
the compressible fluid's viscous stress tensor t is considerably 
more complicated than the formulation - equation (III. 11) - for a 
Navier-Stokes fluid. The structure of equation (III. 9), with its 
characteristic Nabla operators (which have multiple effects on 
the local specific entropy s of the fluid) , contains the 
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asymptotic limiting case 


lim 

<j-> 0 


7 = 0 


£ ^ £c tion- 

behaviour is cloJh, !!!^ '* The demand for sn^h 1 local 
equilibrium On i-ho^i related to Maxwell's nrinri ^ an as ^ Tn P t °tic 
a flow field ^Y® 1 of the Maxweu! B ofi?» lf>le ° £ kinetic 

fS 1 ^S 1 |f. k ^Wc a ^Ubri™ is 
P. 414,. CUrrent densities equal zero (sfe ??ueIdeU , 

Equation fin , . 

SSte ;•« cannot be 

° f the — 


T 


^ + £ s k 3 it (RZ p) v 


(III. 13 ) 

H r ££Vbl€^ fft “ Slt °" dansi .£"!«* q and a convective 

K%%^2n P F ea ^ 

?hf nnea C r S l the 7 «w™ 3 y* n “ *"■ 

proof; even r and e odd ble d TherIn0dynamics and^Sstu^t iVS Part of 
equation for r-r,L° ° rder tensors may not lates wi thout 

Polenicized aqai^rth- 1 ”® 0113 "^ 3 - <=. V Truesde?? ar 5° gether in an 

”!)■ He no? 9 oiry t tracld"?t rlnClple " (*S3Sl 1*1 

Un 1 9, 12 Sfer° Ve 311 its ” at HeTOtTcri d a„bi US Ori9ln ' d “t 
principjf J? a V xaII,ple of the probl^s »n lty - E ? ua ti°n 

follow two sc f al°ars diVe T r h g e e y nCe ° f af£ the 

^SSS. t SrS“i Vith th ^ ^ ev£ - n ouTst°^s e s V \ n en¥ er b 2e - “ h ose 

motion (see Moore, 19 ?!?^.“^ Hevier-stokes^?^ 5 ^ th i 3 
simult^n 2 oisly r oc 2 urina e ? aUSe its " de duction" fail ^ 

constJucUo^Sucras^n 0 - ' ^oSever* 

accepts the known incompatibili ty S of ^his^ f: "onJ 1 

this theoretical fluid with 


14 


the thermodynamics of real materials, one is compelled to 
acknowledge that concepts such as temperature, entropy, heat flow 
density, etc. are not, strictly speaking, compatible. Taking 
this circumstance into consideration and assuming 3 = a for the 
shear viscosity u, it is notable that the Navier-Saint Venant 
motion equation (III. 9) for an incompressible model fluid (v • v 
= 0) contains the Navier-Stokes motion equation as a 
"pathological" limiting case (see Straub, 1988, pp. 141-144). 

Seen from the standpoint of the AT, a notable conclusion can 
be made: for an incompressible fluid with constant shear 
viscosity u, the Navier-Stokes equation of motion, together with 
the incompressibility condition 

7 • v = 0 (III. 14) 

and the compatible initial boundary constraints, represent a 
complete set of equations. The solution - flow velocity v and 
pressure p dependent on location r and time t - is parametrically 
dependent only on the kinematic viscosity v = u/ p = constant 
(and, additionally, on the parameters fixed by the initial and 
boundary constraints) ; to a certain extent individualizes the 
incompressible fluid with a number for comparison. 

In reality things look completely different: for all 

compressible fluids the new quantity 3 (a physically 
corresponding function of u) - equation (III. 10) - is normally a 
function of state, provided the proportionality coefficient t T is 
either a constant or a known function of the density P , the mass 
fraction u k and temperature T. In the latter case this 
characteristic time t T is a property of state like the pressure 
p, the compressibility factor Z and the absolute specific entropy 
s of the polynary single-phase system. 

Even simple examples, however, lead one to expect that this 
coefficient 3 is more likely to be dependent on the class of 
given flows (for example, a shear flow under certain constraints 
such as constant shear and isothermal fluid) than on the material 
properties (see Straub 1988, section 5.2). It is evident that 
this difficulty raises questions about conventional viscosity 
concepts. This is not the place, however to present an answer to 
this problem, particularly since other authors have already made 
some pertinent attempts (see, for example, Nettleton 1987) . If 
necessary, there are no theoretical objections against accepting 
t T as an empirical coefficient for the present. 

Generally speaking, one arrives at an interesting 
conclusion: the Navier-Stokes equation of motion can be 

substantiated only for incompressible model fluids; for real 
compressible materials, one can only define it as addressed by 
Truesdell [9]. 

It is understandable that the Navier-Saint Venant equation 
of motion at present cannot be used for actual projects. 
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Currently . available numerical procedures and computers are still 
too inefficient to handle the work. Therefore it is appropriate 
to search for a first non-trivial approximation of equation 

. The answer naturally depends on which approximations 
one can agree upon. 

The gas kinetics fundamentals of non-uniform and thermally- 
perfect gases lead one to the conclusion that the most important 
prerequisite for a decisive simplification of the problem, the 
inequation 

vT/T « 1 , (III. 15) 

may be valid for the profiles of local temperatures. in equation 
(III. 15) vT relates to the absolute value of the local 
temperature gradients. If one ignores extreme cases, the 
condition (III. 15) is still fulfilled within the linear 
dimensions of the order of magnitude (V m /N A ) 1/3 = 33 a [ 6 ]. This 
means that in most practical cases there are no relevant 
limitations. 

In these cases, both Navier— Saint Venant equation of motion 
(III. 9) and the energy field equation can be simplified further 
in the form: 


c v pDT = - pv • v - 7 • q F + 1/2 

u„ (V • V) 2 

(III. 16) 

pDv = P f - vp v 


(III. 17) 

Pv* P ™ 1/ 2 /X ^ V 

• V 

(III. 18) 

u M : = 

ft 

(III. 19) 


p = RT P 

(III. 20) 


The same can be done with the field equations of the mass 
fractions not shown here. 

Equation (III. 17) agrees formally with the Euler equation of 
motion. Yet the pressure p v by no means corresponds with the 
thermodynamic pressure p according to equation (III. 20). 

Equation (III. 17) describes dissipative processes and can, like 
the Navier-Stokes equation of motion, fulfill no-slip boundary 
conditions. Analogies due to numerical viscosity effects 
obviously are the reason why the Euler equations of motions have 
recently been so successfully applied in industry for simulating 
real flows. 

An equation of motion of the Euler type (and not the Navier 
Stokes equation of motion!) is thus a first step toward 
describing real flow fields of compressible fluids. All field 
properties such as velocity v, specific field force f, mass 
density Pl pressure p, and temperature T are clearly defined by 
the Gibbs Fundamental Equation and related to one another. 
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If one accepts for q F the known proportionality to 7 T 
according to the First law of Fourier 

q ? : = - x 7T = - 7 • [c v u. n f 7T] . (III. 21) 

that is equivalent to assume steady transfer processes (see 
Straub et al. 1987) in connection with the Fourier heat flow 
vector q F . The coupling of heat conductivity \ with the Maxwell 
viscosity (analog to the known gas kinetic relationship x = f c v ^ 
between X and the dynamic viscosity t) of non-uniform gases [6] V 
allows further simplification of equations (III. 16) and (III. 17), 
assuming continued constant material data. If one adds the well- 
known relations 

c V /R = (* - l)' 1 ; A: = p/p = RT (III. 22) 

for perfect gases and ignores the quadratic v • v term in 
equation (III. 16) as well as the field force term in equation 
(III. 17), then these two equations are immediately transformed - 
taking equation (III. 15) into account (see Straub, 1988, p.148) - 
into the forms 


DA = - (* - 1) [A 7 • v - 1/2 v n 1 • 7 A] K\ = c p/c v (III. 23) 

D v = - 7A + 1/2 7(7 • v) ; v n : = u „/ p (III. 24) 

They are of particular interest because they contain only two 
dependent variables - flow velocity v and the scalar potential A, 
which with constant material data can also be easily connected 
with the special enthalpy of a perfect gas - as well as two 
parameters i and To solve them, one needs suitable initial 

and boundary conditions, but not the continuity equations: the 

latter can subsequently be used, when fields A(r, t) and v(r,t) 
are established to calculate the density field p(r,t). 
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IV. 




S”;”“ »• a...i«i h„.,i = 

[®]* Particularly, we shall dpmnnc-i P J d by Chapinan and Enskog 
approximation is equI?aien“to th» » * t ?f t the f Vst order 

fluid flow as well’as gasSJnamics Navler -St°kes equation of 

1V.1. lie Boltzmann Eauation and the Equation o! Transnort 

f(v,r,t), such f that thS l0Cal dlstrlbut i°n function for the gas, 


£ (v, r, t) dv dr 


(IV. 1) 


positiSff/f+^fU Ume a t tiC ThI Joci? g d Ve t° CitY * + at 
a function of the macroscopic properties sinh butlon function is 
mean velocity and temperature P Thus \h P as number density, 
distribution function f(v,r t) at a tS ° f chan ^ e of the 

t by considering a ga? is sub ec^ ?o an P °f ltion r *nd time 

may be a function of ? and t but nS 2 external force mF, which 
a "d t + dt, the velocity V of v ° i 7* Between the time t 
with another, phanges to v -/rd/ and Whl f h d ° es not ?°Hide 

change to r + vdt. There are f (v ? ^H^ P ° Sltl ° n “ vector r will 
time t^ lie-in the volume-element J olecules which at 

range v, dv. After the intlrva^ at if ?h ^ velocities in the 
could be neglected, the same molecules effect of collisions 

compose the set that occudv thp vniV - nd - n ° ot £ ers > would 
veiocities in the range H 15? & “?.! * V * fc ' db ' abd have 

y t iar, dv, the number in this set is: 

f(v + Fdt, r + vdt, t + dt) dvdr . 

equ??I5? 1 ? g r b ? h f5^A?5?i?? k f5?c??o5 e ? d i5 0 0 ??r?' e | olt2 ” ,ann ' s 


df 

-t 

df 

af 

a.f 

+ V < 

1 — + F 

f -i— i m — 

e 

at 

3r 

r > 
(X> 

at 


(IV. 2) 


The quantitiy of -t- defined above is equal to the rate of the 
change, owing to encounters in the velocity-distribution function 
of f at a fixed point, it will appear that !fi is ex pressible as 
an integral involving the unknown function f?* 
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\ 3t / 2 t 

, -f 

with k 12 dk the collision parameters. 

Another important equation may be derived from Boltzmann's 
equation which give the connection between the molecular velocity 
distribution and conservation laws of fluid dynamic equations. 

Let <t> be any molecular property as defined by: 


ff' f' - f f \ k dk dv (IV. 3) 

J [ 1 2 1 2 j 12 2 ' 


0 


1 P 

n i 


?>fdv 


(IV. 4) 


with n being the number density of molecules and mass the density 
given by nm, with m being the mass of a molecule. J 


Multiply Boltzmann's equation by 0dv and integrate 
throughout velocity-space; it is assumed that all the integrals 
obtained $re convergent, and that products such as {if tend to be 
zero and v tends to infinity in any direction (i.e., -©, +©) . 

The result may be written as, 


C3 


i ( 3 f 

3f „ 

dV 

0 | + V 1 

1 — + F 


J lat 

3r 

dv) 


dv = nij( (IV. 5) 


where 


na0 


a.f - 

0 dv 

at 


(IV. 6) 


and the integration limits are from--© to +®. The physical 
significance of ±<p is the measure of the rate of change., in the 
molecular property with velocities in the range of v, dv per unit 
volume at r, t, owing to encounters. 


From Eq. (IV. 5) , we 


' 3f _ 

3 

P 

P 

0 — dv 

= — 

pddv - 

J 3t 

3t t 


J 

f . df 

3 



0V • — 

—i 

1 

0vfdv 

J 3r 

dr 

t 



obtain 

30 3n0 30 

f — dv = - n — 

3t 3 t at 

^ - 3 — 

fv • — dv = — • n0v 

J 3r 3r 


(IV. 7) 


_ 30 

nv - — 
3r 


(IV. 8) 
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dvdw 


(IV. 9) 


1 af 

<P — dv = 
av 


i p 

[*f] 

J t 


U s oo 


U = - 05 


' 30 a^ 

f — dv = -n — - 
J 3v 3v 


Substituting these results into Eq. (IV. 5) we obtain the 
transport equation: 

dn't a — l' a a a ^ 

+ — • njiv - n — + v ♦ — + F • — = n^0 (IV. 10) 

at 3r lat 3r dv) 

This _ equation could be expressed in terms of a peculiar velocity, 
C = v - c 0 , with c 0 being the center of mass velocity. With 
proper mathematical manipulation, we obtain the transport 
equation in terms of the peculiar velocity: 

D 3 d 7 (Dji aji ( Dc 0 ) aji 

— n<j> + n<p — • c 0 + — • n£c - n j — + c — + Jf - 1 • — 

Dt ar ar l Dt dr { Dt ) dc 

df> d ' ) 

c — c 0 = n±<(> (IV. 11) 

dc dr J 

D 

with — being the substantial derivative, which is given by 
Dt 

Da d 

— = — + c 0 • — (IV. 12) 

Dt at ar 


The transport equation (IV. 11) is used for the bridge to 
connect the microscopic theory to the continuum (fluid) theory. 

IV. 2. Conservation Laws of Fluid Flow: Special Forms of the 

Equation of Change of Molecular Properties . 

The conventional form of the governing equations of fluid 
flow could be obtained from the transport equation by assigning 
the proper form of the molecular property <j > . 

d <p D0 d i> 


Let < 

P = 

1, <t>c = 0, — = 0, — =0, — =0, 

II 

O 



ac Dt ar 


Thus 

Eq 

. (IV. 11) becomes 


Dn 


a 


— + 

n ■ 

— * c 0 - 0, 

(IV. 13) 

Dt 


3r 

This 

is 

the equation of continuity. 
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II. Let 0= 

mC ; 

then 

0 = 

0, n0c = 

-4 -4 

ncc = 

IT (i.e. 

pressure 


3 0 


D0 

3 0 


30 


tensor) 

' ~ 

= 0 

r 

= o, — 

= (m) , 

- = o, 

L0 = 0 . 


dr 


Dt 

3c 


3c 


Hence , 

Eq. 

(IV. 

ii) : 

becomes 




Dc 0 

3 


-» 





P 

-i 

• n 

- § F 

r 9 • • 



(IV. 14) 


Dt 


3r 


This momentum equation at the lowest order 
corresponds to the Maxwellian distribution function 
which will lead to the momentum equation for 
inviscid flow. The next higher order will lead to 
the Navier-Stokes equation for both compressible 
and incompressible flow. 

1 _ 3 — 

Case III. Let <f> = E = — mC 2 , then 9 = — nkT, nf>c = q (i.e. 

2 2 


80 D0 

heat flux) , — = 0, — = 0 and, since E depends on 
8r Dt 

—4 

c only through the contribution of the kinetic 

8 0 H 8 0 

energy of translation, — = me; thus — = 0, and 

8c 3c 

8 0 

r -4 -4 -4 — 

n — c = pcc = IT. Since A0 = 0, the transport 
3c 

equation becomes, 


DT 2 


IT 


3 


Dt 3nk l 3r 

this is the energy 



equation. 



(IV. 15) 


IV . 3 . Maxwell ian Velocity Distribution 

Consider a simple gas whose molecules are spherical, 
possesses only energy of translation, and is subject to no 
external forces. If its state is uniform; 
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5 f 
St 


(f'fj - ff ! ) k ! dkdc , 


f ' f ! = ff t Detailed balancing (i.e. at statistical 

equilibrium) 

log f' + log f[ = log f + log 

This shows that log f is a summational invariant for 
encounters. Thus, it must be a linear combination of the three 
summational invariants (conservation laws) : 

log f = I a [i> i£ (i) 

= a (1 > + a (2> -me - a (3 > - me 2 . 


The form of the Maxwellian distribution is 
a ' ’ ‘ e 


( 3 ) 7 

f = „ < o ) 1 / 2 m c 1 


log f = a (1) + m(a (2) u + a t2, v + a (2, w) - 

! V rr 


a 


( 3 ) 


m (u 2 +v 2 +w 2 ) 


= loga (0) - a 1 


3 ) 


( 2 > 


a. 


u - 


( 3 ) 


a 


( 2 ) 


V - 


„ ( 3 ) . 
a j 


( 2 ) 


W - 


a„ 


„ ( 3 ) , 
a ) 


a 


( 2 ) 


a 


( 3 


By definition 


n = 


fdc = a 


_ „ { o ) 


-a { 3 J 1 / 2 me 


2 ' 


- „ < o ) 

— a 

1 c 2' e -a t3) 1/2 mc 2 ' dcV 

1 P 

sin 6d9 

i 


L 


d^ 


= a (0) 


r 2 7T 


\ 3 / 2 


\ ma ( 3 1 } 
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n 


nc 0 = 


cfdc 


a 


( 3 ) 


+ c 1 


Va 


f dc ' 


' ( 2 > 
a 


= n 


a 


( 3 ) 


+ a 


( 0 ) 


{ 3 ) 


. 2 , 


a l/2mc' c x dc 


a ( 2> 


c o = 


a 


( 3 ) 


f = a (0> e _a ' 3 * 1 / 2mc 2 


( ma 


= n 


( 3 ) N 3 / 2 


2 TV 


( 3 ) 

1 / 2 m c 


(IV. 16) 


Using, 


— kT = — me = 
2 2 

m 

2 


m 


2n 


'ma ( 3 1 1 

p 

l 2 TT , 



c 2 fdc 


C 2 e -a< 3 ’ 1 / 2mc ' 2 d< ; 


we obtained a 


( 3 ) _ 


2a 

1 

kT 


( 3 ) 


, thus Eq (IV. 16) becomes 


f = n 


( m 


12 7ikT/ 


3/2 


2 k T 


This is the Maxwellian distribution. It is well-known 
[10,11] that gas has Maxwellian distribution, it corresponds to a 
non— dissipative state which is equivalent to the Euler equation 
of motion of fluid dynamics. 


IV. 4 . Remarks 


It should be noted that this formalism for the derivation of 
conservation laws is only valid for the dilute gas, because of 
the limited treatment of the Boltzmann collisional integral, in 
which only the binary collisions are considered. However, this 
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formalism does not give us a clear track to chase these physical 
processes as expressed in kinetic theory. For example, the gas 
behavior follows a Maxwellian distribution; with it lowest order 
it leads to the Euler equation of motion, whereas the next higher 
order behavior of the gas leads to the Navier-Stokes equation of 
motion. A summary of these comments is depicted in the Table 
below. 
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V. PRIGOGINE'S NON-EQUILIBRIUM THERMODYNAMICS: 
CONSERVATION LAWS* 

It has been understood that the description of a "mechanical 
system" is in terms of the coordinates and momenta of the 
molecules, or in terms of its wave function. However, such a 
description, when applied to systems of interests in chemical 
physics and fluid dynamics such as combustion chamber flow, leads 
to great practical and conceptual difficulties, because it is 
impossible to say 10 2 3 molecules in a macroscopic system. In 
order to remedy these difficulties, the combined thermodynamic 
and fluid dynamic methods are used to provide us with a "reduced 
description," as pointed out by Prigogine, "simplified language" 
with which to describe macroscopic systems. In this chapter, we 
shall briefly summarize Prigogine 's concept through the 
derivation of conservation laws and balance equations which are 
equivalent to the equation of motion of combustion chamber flows 
being used for current CFD simulation. 


Let us consider a system of volume V limited by the surface SI. 
We wish to follow the time change of the intergral 


I(t) 


fdV 


(V.l) 


This intergral is extended over the volume V of the system we are 
interested in. The surface fl is assumed to be at rest. In the 
terminology usually adopted in thermodynamics, I(t) is an 
extensive quantity. For instance, it may be the mass or the 

energy of the system. On the contrary, f(x, y, z, t) is an 

intensive quantity, which does not depend on the system as a 

whole. It corresponds to the volume density associated to I and 

may be represented by the functional derivative 

51 

f E — (V. 2 ) 


The change can be expressed by 

31 

— = P[I] + (6 [I] 

at ( V . 3 ) 


*The materials presented in this chapter are based on a text book 
entitled Thermodynami c Theory o f Structure, Stability and 
Fluctuat ions , by P . Glansdorff and I. Prigogine: Wiley- 

Interscience, 1971. 
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The first term on the r.h.s. corresponds to the production per 
unit time of the quanity I inside the volume V. It can be 
written as a volume intergral 


P[I] 


n 


cr [ I ] dV 


( V. 4 ) 


where <r [ I ] denotes the source of I per unit time and unit volume. 


The second term in the r.h.s. of (V.3) represents the flow of the 
quantity I through the boundary surface ft. It can be written as 
a surface intergral 


0[I] 



[I]dft 


( V. 5 ) 


This formula introduces the density of flow j [ I ] associcated with 
I; j n [I] is its projection along the inside normal to the 
surface. For simplicity, we shall also use the term flow instead 
of density flow (or current, or flux). From equations (V.4) and 
( V. 5 ) we obtain the so-called balance equation corresponding to 
the extensive variable I as: 


ai 

at 


1 

o' [ I ] dV + 

l 


n 


j n C i ] an 


(V.6) 


It is sometimes useful to write the balance equation (V.3) or 
(V.6) in the symbolic form 

dl = d | I + d e I ( V . 7 ) 

where djl represents the source term and d e I the flow term. One 
may write as well 


d i I = dl + (~d e I) (V. 8 ) 

In this form, we see that the source djl contributes on one side 
to the time change of I and on the other to the flow from the 
system to the outside world (-d e I). It must be emphasized that 
only dl is in general a total differential of the state 
variables . 


Equality (V.6) has to be valid whatever the volume V. Therefore, 
the application of Green's formula gives us directly the balance 
equation in local form 

df 

— = <r [ I ] - div j [I] (V. 9 ) 

at 
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One of the advantages of this formulation is that all 
conservation equations can be expressed by the statement that the 
source term corresponding to a conserved quanity vanishes. For 
example, if I represents the total mass system on has 

I = M 

f becomes the mass density p 

f = P 

and conservation of mass is expressed by the relation 

<r[M] = 0 (V.10) 

The time change of the density p is then, apart from the sign, 
equal to the divergence of the mass flow. Moreover, the mass 
flow is clearly 

j [M] = pv (V. 11) 

where v is the velocity of matter. Using (V.10) and (V.ll), we 
therefore obtain the classical continuity equation 

dp 

— + div pv = 0 (V. 12 ) 

at 

Likewise, conservation of total energy U (first law of 
thermodynamics), and of total momentum Q (in the absence of 
external forces) , may be expressed as 

<r[U] = 0; <r[Q] = 0 (V. 13 ) 

The source term for entropy S plays a special role as the second 
principle of thermodydnamics and postulates the inequality: 

<t[S] > 0 -- ( V. 14 ) 

Entropy is not a conserved quanity, but increases as the 
result of irreversible processes included in the source term. It 
is only for reversible processes that the change of entropy is 
entirely due to entropy exchanges with the outside world. 

Let us go back to the density of flow j[I] in the balance 
equation (V.9). In general, we are concerned not only with a 
convection flow such as (V.ll) but also with a conduction flow 
j cond (say), which occurs even in a system at rest. We have 
therefore 

j — j c o n d 1 c onv ~ Icond f V ( V . 1 5 ) 
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For example, the heat flow W is the conduction current 
associated with internal energy. We shall consider other 
examples later on in this chapter. 

In many cases, it is useful to introduce the so-called 
substantial or hydrodynamic derivative 


d3 a 

— = — + I Vi (V. 16) 

dt at j dx i 


The continuity equation (1.12) leads then to the equalities: 


d p 3 3 

P — = p — + E pv j 

dt at j 3Xj 


3 ( pi>) 3 

+ l ( pV ; P ) 

at j ax. 


( V. 17 ) 


This relation may be applied to an arbitrary intensive variable 
o(x,y,z,t) = <p (x x ,x 2 ,x 3 , t) . 

A summary of these results are shown in Table IV. Table IV 
also provides a comparison with solutions from Straub's AT 
method. 

The details for the derivation of conservation laws can be 
found in the book entitled "Thermodynamic Theory of Structure, 
Stability and Fluctuation" by Glansdorff and Prigogine (1971) [7]. 
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VI. TRUESDELL' S RATIONAL THERMODYNAMICS: A CLASSICAL REPROACH 


The approach taken by rational thermodynamics researchers is 
to seek concepts and principles sufficient to describe and 
calculate temperature fields in bodies subject to deformation and 
heating. Rational thermodynamics has pointed its attention at a 
process occurina in time . as opposed to the classical approach of 
Gibbs who has sought in statistical mechanics and from it has 
constructed a rational foundation for the statics of hotness, 
heat, and equilibrated forces. Recent applications of rational 
thermodynamics were the propagation of stress waves, the effects 
of severe deformation upon a body's ability to conduct heat and 
of heating, upon a body's response to stress and strain, and the 
nature of energetic and diffusive transfer, incident upon change 
of shape. Because such effects notoriously differ in bodies of 
different materials, an adequate theory must distinguish 
materials by constitutive properties, openly and explicitly 
stated. The constitutive function of a body must make it 
impossible for that body ever to violate the fundamental axioms 
of thermodynamics. The fundamental theory of rational 
thermodynamics is to reproach that the concepts and assumptions 
of the pioneers such as Kelvin, Clausius and Gibbs, selected and 
refined and extended, can be made the basis of a treatment that 
reflects simply the physics of heat and work and is not only 
rigorous, but also broad enough to cover irreversible processes 
in deformable continua. 

The treatment of dissipative phenomena highlights the spirit of 
rational thermodynamics of Truesdell. In terms of describing 
fluid flows, rational thermodynamics rests upon some prior 
concepts of mechanics which may mean the ordinary theory of a 
deformable body. The concepts of mechanics were further extended 
to allow for diffusion and chemical reactions as well. 

In developing a thermodynamic theory for deformable bodies, 
Truesdell starts with the classical theory of continua. Let a 
body B occupy a domain of space at each time t. These domains 
are its configurations. The mapping of one such domain or 
reference configuration K onto the present configuration is 
called the motion of B with respect to that configuration: 

X = X K (X, t) , - CO < t < CD. (VI. 1) 

Here x K is the mapping that carries the place X in the reference 
configuration K into the place x at time t. The places X 
identify the particles; thus we may speak of "the particle X" 
without fear of confusion. The velocity x is the rate of change 
of position of a particle: 

x = a t x K (X,t) = x ( x , t ) . (VI. 2) 

A mass, which is a non-negative measure M, is defined once and 
for all over measurable subsets of the body and is also assumed 
to be an absolutely continuous function of volume. Hence a non- 
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negative mass-density p ( • ) exists almost everywhere at each time, 
and at points where it and the velocity field x( • ) are 
sufficiently smooth 

p + p div x = 0 . (VI . 3 ) 

The linear momentum m and rotational momentum fl xo of B in its 
present configuration are defined by 


m e 


x dM, 


fl 


X 0 


( x - x 0 ) A x dM , 

“ B 


(VI. 4) 


where x 0 is a fixed place and where the integration is based on 
mass. Euler's laws of mechanics assert that for every body the 
time rates of change m and li xo of these quantities are equal to 
the total force f and total torque F xo acting upon that body in 
its present configuration: 

m = f, ii X o = F x0- (VI. 5) 


In the simplest kind of continuum mechanics, the force f is 
assumed to be the sum of two forces of special kinds: the total 
body force f b and the total contact force f c . The total body 
force is an absolutely continuous function of mass, while the 
total contact force is an absolutely continuous function of 
surface area: 


f = fc = f b / 

i p 

= tdA + bdM. (VI. 6) 

J 3B J B 

The field b, defined in the interior of B, is the body force per 
unit mass; it represents forces acting at a distance, e.g., those 
of gravitation and electromagnetism. The traction t represents 
the action of neighboring parts of material upon one another, or 
any forces applied in any way to the boundary 3B of the body B in 
its present configuration. According to a theorem of Cauchy, the 
traction is delivered by the stress tensor T: 

t = Tn, (VI. 7) 

n being the outward unit normal to 3B. By combining (VI. 4), 

(VI. 5), (VI. 6), and (VI. 7), on the assumption that the various 
fields occurring in those equations must be sufficiently smooth , 
it is easy to derive a local equivalent to the principle of 
linear momentum for continua: 
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div T + P b = px, (VI. 8) 

called Cauchy's first law of motion. 

If we restrict attention to the case when all torques are 
moments of forces, and so 


X 0 


8B 


(x - x 0 ) A t dA + 


(x - 
9B 


X 0 , A b dM, (VI. 9) 


then another theorem of Cauchy asserts that the stress tensor T 
is symmetric: 


T = T t (VI. 10) 

This is Cauchy's second law of motion. 


The kinetic energy K of B in the present configuration is 
defined by 


K 


1 
2 tl 


x 2 dM , 


(VI. 11) 


while the power P is the rate of working of all the forces acting 
on B: 


P 


t 

J 3B 


x dA + 



x dM. 


(VI. 12) 


Now, according to a simple theorem -first proved by Stokes, the 
net working W at a time T on the body B is given by the 
difference of P, the mechanical power, and K, the kinetic energy, 
thus: 


W = P - K, 


w 

J B 


dV, 


where 


w = T • grad x, 


(VI. 13) 


(VI. 14) 


the dot indicating the inner product in the vector space of 
second-order tensors (X • Y = tr (XY T )). The scalar w is the net 
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working per unit volume; it is often called the stress power. 

These few formulae provide a specfic mechanical framework upon 
which a structure of thermodynamics may be raised. Here, of 
course, all quantities are given in terms of time-dependent 
spatial fields, not merely functions of time alone, and so we 
cannot blindly apply the results of Lecture 1, but we may proceed 
in a parallel way. A number of things I shall say will be valid 
equally in more general kinds of mechanics, where a result of the 
form (VI. 13) holds with a density w which is given by an 
expression more complicated than (VI. 14), allowing, for example, 
for the action of couple stresses, but, so as to have a fully 
definite and explicit mechanical background. 

The first law of thermodynamics, expressing the balance of 
energy is [9]: 


E = W + Q 

where E is the internal energy of a body, W is the net working on 
it, and Q is its heating. This axiom we now set alongside 
Euler's laws (VI. 5) as basic for all theories relating forces and 
energies to motions. Just as we have specialized Euler's laws by 
assuming (VI. 6) so as to obtain forms appropriate to continuum 
mechanics, so also we may specialize the principle of balance of 
energy in a corresponding way. First of all, we take the net 
working W to be that calculated according to continuum mechanics 
and hence given by (VI. 13) as an absolutely continuous, additive 
set function with density (VI. 14). Parallel to the assumptions 
(VI. 6) and (VI. 9) regarding the resultant force and resultant 
torque, we lay down the assumption that the heating Q is the sum 
of the two heatings of special kinds: the body heating Q b and 

the contact heating Q c ; the former being an absolutely continuous 
funcntion of mass, and the latter, an absolutely continuous 
function of surface area: 

Q = Q c + Q b / 

1 p 

= q dA + s dM. (VI. 15) 

t 3 B i B 


The volume density s is called the heating supply; the surface 
density q, the influx of heating. The two kinds of heating, Q c 
and Q b , are often said to describe conduction and radiation, 
respectively, but no such specific connotation need be made in 
general. Finally, the internal energy is assumed to be an 
absolutely continuous, additive set function: 


E 


e dM, 

J B 


(VI. 16) 
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the density e being called the specific internal energy or, for 
short, the energetic. 

Under the assumption (VI. 16), from [9] and (VI. 15) Stokes 

derived by Cauchy's method a result inferred earlier in a major 
special case by Fourier, namely, the influx of heating is 
delivered by the heating-flux vector h: 

q = h • n, (VI. 17) 

n being the outward unit normal to 3B. For sufficiently smooth 
fields, by (VI. 17), (VI. 16), (VI. 15), (VI. 14), and [9] it is easy 
to show by use of the divergence theorem and (VI. 10) and (VI. 8) 
that 


pf = w + div h + pS. (VI. 18) 

This differential equation is sometimes regarded as a local 
statement of the "first law of thermodynamics", although in fact 
it is a consequence of a rather subtle sequence of field 
definitions along with Cauchy's laws of motion (VI. 8) and (VI. 10) 
as well as the general balance of energy [9]. Major special 
cases of this differential equation were derived, under various 
special hypotheses, by Fourier, Kirchhoff, and C. Neumann. 

In the thermodynamics of homogeneous processes the axiom of 
irreversibility takes the form of the Clausius-Planck inequality 
[9]: 


6H > Q, 

where 9 is the absolute temperature and H is the calory. We must 
now introduce the concepts of temperature and calory in continuum 
mechanics. The former offers no difficulty. The temperature 9 
is given by a positive-valued temperature field, 

9 = 9 (x, t) _> 0, (VI. 19) 


defined over the persent configuration of the body B. As in 
statistical mechanics, also in continuum mechanics the reciprocal 
of the temperature, which may be called the coldness J is often 
more convenient for the mathematical theory: 

1 

J E -. (VI. 20) 

9 


It is natural to assume that the calory, like the internal 
energy, is an absolutely continuous function of mass: 
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H 


1/ 


(VI. 21) 


dM 


The density r, , which is usually called the specific entropy, I 
prefer to name the coloric. This name has the virtue of dragging 
a red herring across the path of those who lay down police-edicts 
concerning when and how "entropy" is to be "defined" . As they 
back off in horror at revival of the long-dead caloric theory of 
heat, we can return to our concrete, mathematical theory, safe 
behind our two negations of authority. 

Caloric is introduced to represent gross dissipation. While 
coloric may flow from place to place and thus increase or 
decrease locally, its total, the calory H, is subjected to a rule 
favoring increase. We may regard this vague statement as 
expressing the content of Clausius' assertion: 

H > 0 when Q = 0. (VI. 22) 

The pioneers of thermodynamics regarded this assertion as a 
general axiom, not restricted to any special type of body or 
process, but they were unable to implement it by mathematics such 
as to give it a concrete form. Their various "second laws" sound 
more like warnings or threats than principles of a rational 
science. In modern continuum mechanics, because (VI. 15) holds, 
we may express Clausius' assertion thus: In a body B where 

conditions are so balanced that 


q dA 

J 9B 


+ 


i 

s dM = 0 , 

J B 


(VI. 23) 


the calory cannot decrease: 


•n dM > 0. (VI. 24) 

J B 


(A body B is said to be in a "calorically isolated" or 
"insulated" configuration if q = 0 on 9B while s = 0 in the 
interior of B. For such a body, of course, (VI. 23) holds.) 

While this statement expresses some idea about irreversibility, 
it does not lead to a unique generalization of the Clausius- 
Planck inequality [9] to deformable bodies. In particular, since 
temperature does not appear in it, it affords no connection 
between the ideas of heat and temperature. In this lecture I 
shall motivate the principle now regarded by most of the 
presently productive theorists, though not by all, as being a 
correct dissipation inequality. To this end I shall present some 
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simple aspects of dissipation in familiar examples. 

The internal dissipation 5 is defined as being the amount by 
which the increase of caloric multiplied by the temperature 
exceeds the local heating: 


1 

5 = 9n - ~ (div h + P s) . (VI. 25) 

P 


While 5 is defined in terms of the quantities associated with 
temperature and heat alone, the equation of energy balance 
(VI. 18) permits us to interpret 5 also as the amount by which the 
increase of caloric times temperature exceeds the increase of 
energetic not produced by working: 


5 = e-n - I 6 - - . (VI. 26) 

l p ) 


If we were to transfer to the local description the ideas that 
motivate the Clausius-Planck inequality [9], we should require 
that this increase be non-negative: 

5 > 0. (VI. 27) 


We shall name (VI. 27) the Planck inequality. 4 For reasons to 
appear later, we shall not impose it as an axiom. Instead, in 
the remainder of this report, we shall first calculate 5 for the 
most familiar kinds of materials and then note some consequences 
that follow when indeed 5 > 0. 

First, for the ideal gas of the Euler-Hadamard theory, 

T = -pi (VI. 28) 


and also [ 9 ] , 


e = « (v, i)) , P = -3 V f , 0 = 3 , 

where v = 1/p. Hence by (VI. 14) and (VI. 3) 

w p 

- = - - tr (grad x) = -pv, (VI. 29) 

P P 

and so by (VI. 26) 
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5 = (3,0 1 + (3 V 0v 
= 0 . 


(VI. 30) 


Thus the internal dissipation of the ideal gas is always null. 

The same result holds for the general thermo-elastic material, 
which includes as a special case the ideal gas. 

Consider now the linearly viscous fluid of the Navier-Stokes 
theory. For it, a viscous stress V is superposed on the Eulerian 
pressure: 


T = -pi + V. (VI. 31) 

and this V is assumed to be a homogeneous, linear, isotropic 
function of the stretching D: 


2D E grad x + (grad x) T , 


(VI. 32) 


and 


V = X (trD) 1 + 2 uD, (VI. 33) 

the pressure p still being related to the internal energy by [9]. 
Hence by (VI. 26) 


p 5 = V • D, 

= X (trD) 2 + 2u tr D 2 . (VI. 34) 

Thus the internal dissipation of the Navier-Stokes fluid is the 
working of the viscous stress V. The requirement 5 > 0 is 
standard in this theory. Moreover, it is imposed identically in 
D. That is, in the theory of viscous fluids we demand that the 
viscous stress shall never give out work in any motion we might 
imagine 5 , but it may use up work, and generally it will. It is 
easy to show from (VI. 34) that the condition 5 > 0 for all D is 
equivalent to the classical inequalities of Duhem and Stokes: 


u > 0 3 X + 2 u. > 0 . 


(VI. 35) 


The concept, described so far in Lecture 2 of [9] by Truesdell, 
contains the important aspects of rational thermodynamics and its 
application to fluid flows. 
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VII. Discussions and Concluding Remarks 


In this study, we have surveyed a number of theories in 
connection to the description of combustion chamber flows. 
Particular attention was paid to the Falk-Ruppel theory of 
'Gibbsian thermodynamics' and Straub's 'Alternative Theory' (AT) 
in addition to Prigogine's theory of 'Non-eqilibrium 
Thermodynamics', Truesdell's 'Rational Mechanics', and Boltz- 
mann's 'Kinetic Equation'. Whenever possible, comparisons were 
made with the conventional Navier-Stokes equations and Straub's 
AT. Table I shows side by side the equations of motion from 
Navier-Stokes and Straub. Table II compares respective 
formulations of Falk-Ruppel and Straub. Wu's formalism, deduced 
from Boltzmann's kinetic equation, is compared with the classical 
Navier-Stokes equation in Table III. Finally, the comparison 
between Prigogine's non-equilibrium thermodynamics and Straub's 
AT is shown in Table IV. 

A direct correspondence between Falk-Ruppel ' s theory and the 
Navier-Stokes equation could not be established at this time. 
However, Boltzmann's kinetic equation with certain constraints, 
extended by Wu's formalism, permitted the derivation of the 
classical Navier-Stokes equation. 

During this investigation, no conclusive results could be 
established. However, some comments about the application of the 
classical Navier-Stokes equation to combustion chamber flows can 
be stated. Before addressing these comments, some important 
characteristics, which are fundamental in combustion chamber 
flows, are summarized below: 

(i). Many species exist in the spray/combustion domain 
under wide ranging pressure and temperature conditions. These 
species could be present also in different phases. This 
identifies the combustion chamber flow as a multi-flow and multi- 
phase system. 

(ii) . During the combustion process the multi-species, 
multi-phase flows go through states of non-equilibrium 
thermodynamics, i.e. each species exhibits its own number 
density, temperature, and velocity. 

( iii) . Since the flow conditions follow non-equilibrium 
thermodynamic states, the associated chemical reactions will also 
observe non-equilibrium behavior. 

With these specific flow characterisics in the combustion 
chamber, it is obvious that the classical Navier-Stokes equation 
becomes theoretically inadequate to describe the complex flow 
system. We realize that the modified Navier-Stokes equations are 
being used currently for combustion chamber flow simulations. 

But these modified equations have been derived on an 'ad hoc' 
basis. To illustrate the inadequacy of the present Navier-Stokes 
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equation for combustion chamber flow modeling, we offer the 
following few comments without exhaustion: 

(i) . It is known that the Navier-Stokes (NS) equation is 
derived for a single-fluid flow. Although the NS equation has 
been extended to multi-fluids with multi-phases, the equation of 
state for ideal gases is used to describe this fluid system, 
which we all know is not the case. 

(ii) . As demonstrated by Wu [10, 11], the radiative 
gasdynamics flow under non-local thermodynamic equilibrium (non- 
LTE) conditions has led to a number of additional transport 
coeficients. These additional transport coefficients will have 
an effect on the momentum and energy transport. Such phenomena 
should be expected for combustion chamber flows. 

(iii). If the momentum and energy transport were modified, 
the flow field will be changed also. Thus, it will affect the 
design and performance parameters of the combustion chamber. 

From these remarks, it is obvious that the governing equations 
(usually called NS equations) deserve a great deal of attention 
to better understand and model the physics of the combustion 
chamber flow. Therefore, we conclude with the following 
recommendations: 

(i). Investigate the potential connection between Falk- 
Ruppel ' s theory of Gibbsian thermodynamics and the Navier-Stokes 
equations . 

(ii). Examine whether Straub's claim from his AT could 
contribute to the improvement of the modeling of combustion 
chamber flows. 

(iii). Conduct a rigorous investigation for a conclusive 
comparison of the NS equations with the respective governing 
equations from the theories developed by Faulk-Ruppel , Straub, 
Prigogine, Truesdell and Boltzmann _as extended and demonstrated 
by Wu [10, 11] . The outcome will establish a reliable foundation 
for the modeling of combustion chamber flows. 
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Table I. 


Comparison Between Classical N-S and Straub's AT 


N-S 


AT 


Assuming start point is Assuming start point is Newton's 
Newton's law law 

dp 

dt 


p = mv -v momentum 


F dv 
v dt 


-=-^P+7*T+pf 


dv _ __ 

p = - 7p + 7 • T + pf 

dt 


No change 


dp 

dt 


p = mv + mp 


F dv d p 

v dt dt 


— = — 7p + 7 • T + or 
v 


dv _ _ dp 

P = - 7 p + 7 • T + pf - f> 


dt 


dt 


dp d p4> _ _ 

P — = + 7 ♦ p ? v 

dt at 


7 

No change 

dv _ _ 3 P <f> 

p = — 70+7 • T + c f — — 7 * p 0 V 

dt at 



assuming P = 1/2 (1 - v) , 1: impulse 

8 

No change 

Dv _ _ 1 9 - 

p = -7p+7‘T+ff + - (pV-pl) 

Dt 2 at 



1 

+ 7 ' p — (V— i) V 

2 


result 

result 

9 

DV _ 

p = — 7p + 7*1 + p f 

Dt 

Dv ___ __ 1 a _ 

p = -7p+7'T s + pf T ( pV— pi) 

Dt 2 at 



here t s = r+— pv (v — 1) 
2 



Remarks: 

10 


1. What does t s mean? 



2. How does one_derive the 
expression r s in AT? 
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Table II. Comparison Between Falk-Ruppel and Straub 


Falk-Ruppel Theory 

Straub 

1 

Gibbsian Thermodynamics 

Alternative Theory (AT) 

O 

Momentum 

Momentum 

Zlr 

p = m v 

p = mv + mp (r, s, p, w k ) 

3 

No difference with 
Newton's Equation of 
Motion under certain 
assumptions . 

There are differences with 
Newton's Equation of Motion 

3 

i.e. an extra term — ( p<f>) 

3 1 

4 


He got his "AT" from Falk's 
dynamics. 

5 

Entropy production is 
the central problem 

He dealt with AT using the 
entropy production 

6 

Reference 
G. Falk-W. Ruppel 
Energy and Entropy (1976) 
Mechanics, Relativity, 
Gravitation (1977) 

Reference 
D. Straub (1989) 

Thermof luiddynamics of Optimized 
Rocket Propulsions 
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Table III-B. Comparison Between Wu's Formalism and 

Classical Navier-Stokes Equation - Continuity 


Wu's Formalism [10, 11] 

Classical Navier-Stokes 

Wu's Formalism on the Basis of 

Classical N-S Equation 

the Boltzmann Equation 


Continuity 

Continuity 

Equation (6.4) of Wu 


dp d 

dp a 

— + — . (pC 0 ) = 0 

— + — • ( pu) = 0 . 

at ar 

at dr 

i 


Eq. (2.22) C 0 = - E pj Uj 


P j 


(2.24a) Uj — C 0 + 3 j 


If all the species have the same 


velocity, namely at LTE. 


=> =0 


Uj = c 0 


pC 0 = Ij P] iij = ( Epj ) u = pu 


d P a 

dp d 

+ • (pU) = 0. 

— + — ■ ( pu) = 0 . 

at ar 

at dp 
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Table IV. Comparison Between Prigogine' s Non-Equilibrium 

Thermodynamics and Straub's AT 


PRIGOGINE 

STRAUB 

1 

Dv 

P = pf - 7 (PS - T) 

Dt 

Dv _ . 1 3 

P = p f - 7 (p 5 T ) + ( pV— pi) 

Dt 2 at 

2 

Micro theory (non- 
equilibrium statistical) 

(A) Starting at Liouville 
Equation 

(B) Parallel to quantum 
Mechanical method 

AT is based on Prigogine 's micro 
theory, but he referred only 
to Prigogine 's book "Chemical 
Thermodynamics" in which the 
micro-theory was not involved. 

3 

What remains unclear in 
Boltzmann's derivation 
is the range of validity 
of his equation. 

The macroscopic laws must be 
throroughly revised according 
to Prigogine' s results. 
Justifying them with the 
Boltzmann equation is untenable. 

4 

In chemical thermody- 
namics, he described the 
entropy production. 

He used the entropy products 
in his AT. 


References : 

Reference : 

5 

1. Prigogine (1962) 
Non-equilibrium 
Statistical Mechanics 

D~ Straub (1989) 

Thermof luiddynamics of 
Optimized Rocket Propulsions 


2. Prigogine (1964) 

Chemical Thermodynamics 



3. Glansdorff and 
Prigogine (1971) 
Thermodynamic Theory of 
Structure, Stability 
and Fluxation 
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APPENDIX I 


Meeting Minutes 
Gibbsian Thermodynamics 

On January 9 to 14, 1991 Professor Wolfgang Ruppel from the 
University in Karlsruhe/Germany visited UAH and MSFC to present 
and discuss a new concept he defines as "Gibbsian 
Thermodynamics." This concept has been proposed in two books 
with Prof. Falk and Prof. Ruppel as authors, written in German 
and entitled Energy and Entropy and Mechanics . Relativity. 
Gravity. 

Professor Ruppel familiarized us with the novel Gibbsian 
Concept, and he pointed out the simplicity of the new concept as 
well as many advantage associated with this approach. 

The concept was accepted in its abstract formulation, but 
questions with respect to some applications were not resolved. 
Prof. Ruppel will give us his answer shortly, after he has looked 
at the problems from the different points of view presented at 
the meetings. 

In general, the new concept appears to be superior to the 
classical thermodynamics approach. The presentations and 
discussions have only covered a limited number of topics due to 
the limited time available. It also takes time to comprehend 
this new approach. However, a dialog has been initiated and 
should be continued. A translation of the two books in question 
or a new report, which is oriented to the evaluation of existing 
equations of motion, will be of advantage. 
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APPENDIX II 


Workshop Report 

Thermo-Fluid-Dynamics Processes in Rocket Thrust Chambers 

The workshop was held at the Mechanical Engineering 
Department of The University of Alabama in Huntsville from 
November 11 - 15, 1991 in Engineering Building Room 102. The 
principal investigator for the MSFC contract is Prof. S. T. Wu. 

This workshop addressed two topics: (a) thermodynamic 

equilibrium calculation for rocket engines, and (b) critique of 
the Navier-Stokes equation, which are discussed in reference (6) 
by Prof. Straub. 

(a) Both, the NASA ODE program and the Munich Method code 
carry now an option to consider finite area combustion effects. 
The solution process in both approaches is based on the 
identified Lagrange function, and the predictions are close. The 
claimed preference of the maximization of entropy over the 
minimization of energy was not accepted. Although comments in 
Ref. (6) relate to an earlier ODE program, documented in NASA SP 
273, there are other specific harsh unjustified claims and 
erroneous results from an inadequate earlier Munich method 
program. To clear the ODE program from these allegations, the 
appropriate steps shall be initiated by Prof. Straub to correct 
this situation. 

(b) Proof of a new 'Navier-Saint Venant' equation validity, 
resulting from Prof. Straub's proposed Alternative theory, was 
not presented and, therefore, the equation was not accepted as a 
preferred scheme over the Navier-Stokes equation at this time. 
Some interesting results, however, should be surveyed. A 
document, showing the derivation of the novel equation, should be 
prepared in English to assure a -tent ion and relevant studies. 
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